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Abstract 

We study the eigenvalue behaviour of large complex correlated Wishart matrices near 
an interior point of the limiting spectrum where the density vanishes (cusp point), and 
refine the existing results at the hard edge as well. More precisely, under mild assumptions 
for the population covariance matrix, we show that the limiting density vanishes at generic 
cusp points like a cube root, and that the local eigenvalue behaviour is described by means 
of the Pearcey kernel if an extra decay assumption is satisfied. As for the hard edge, we 
show that the density blows up like an inverse square root at the origin. Moreover, we 
provide an explicit formula for the \/N correction term for the fluctuation of the smallest 
random eigenvalue. 
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1 Introduction 

Empirical covariance matrices are natural random matrix models in applied mathematics 
and their study goes back at least to the work of Wishart [35]. In the large dimensional 
regime, where both the size of the observations and of the sample go to infinity at the same 
speed, Marcenko and Pastur provided in the seminal paper [25] the first description of the 
limiting spectral distribution for such matrices, see also [31]. Eor instance, this limiting 
distribution has a continuous density on (0, oo); its support is compact if the spectral norm of 
the population covariance matrix is bounded; it may include the origin and may also present 
several connected components. 

Afterwards, attention turned to the local behaviour of the random eigenvalues near points 
of interest in the limiting spectrum, like positive endpoints (soft edges), see e.g. [20, 22, 3, 14, 
17], interior points were the density vanishes (cusp points) [26], or the origin when it belongs to 
the spectrum (hard edge) [15, 17]. Complex correlated Wishart matrices, namely covariance 
matrices with complex Gaussian entries, play a particular role in such investigations since 
their random eigenvalues form a determinantal point process. Indeed, for determinantal point 
processes, a local asymptotic analysis can often be performed by using tools from complex 
analysis such as saddle point analysis or Riemann-Hilbert techniques. In the more general 
setting of non-necessarily Gaussian entries, one then typically shows that the local behaviours 
are the same as in the Gaussian case by comparison or interpolation methods, see e.g. [23, 24]. 

Eor complex correlated Wishart matrices, a fairly complete picture of the local fluctuations 
at every edges of the limiting spectrum has been obtained in the recent work [17], provided that 
a regularity condition is satisfied. This condition essentially warrants the local fluctuations to 
follow the usual laws from random matrix theory. Eor instance, if one considers a soft edge of 
the limiting spectrum, then this regularity condition ensures that the limiting density vanishes 
like a square root at the edge, and the fluctuations of the associated extremal eigenvalues follow 
the Tracy-Widom law involving the Airy kernel. As for the hard edge, when it is present, the 
fluctuations are described instead by means of the Bessel kernel. 

The aim of this work is twofold. Eirst, we investigate the local behaviour of the eigenvalues 
near a cusp point which satisfies the regularity condition: We show that the limiting density 
vanishes like a cube root near the cusp point (hence justifying the name) and, under an 
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extra assumption on the decay of a speed parameter, we establish that the eigenvalues local 
fluctuations near the cusp point are described by means of the Pearcey kernel. 

Our second contribution is to strengthen the results of [17] concerning the local analysis at 
the hard edge: We show that the density behaves like an inverse square root near the origin, 
and we provide an explicit formula for the next-order correction term for the fluctuations. 
This last result is motivated by the recent work [12] by Edelman, Guionnet and Peche where 
they conjecture a precise formula for the next-order term for the non-correlated Wishart 
matrix, a conjecture then proven right by Bornemann [6] and Perret and Schehr [29], with 
different strategies. Our result hence extends this formula, with an alternative proof, to the 
more general setting of correlated Wishart matrices. 

The reader interested in a pedagogical overview on the results from [17] and the present 
work may have a look at the survey [18]; it also contains further information on the matrix 
model and lists some open problems. 

Let us also stress that, at the technical level, the study of this matrix model shares similar 
features with the study of the additive perturbation of a GUE random matrix [9] , and random 
Gelfand-Tsetlin patterns [10, 11], although each model ultimately brings up its own share of 
technicalities. 

We provide precise statements for our results in Section 2, and then prove the results on 
the density behaviour in Section 3, the cusp point fluctuations in Section 4, and the expansion 
at the hard edge in Section 5. 

Acknowledgements. The authors are pleased to thank Eolkmar Bornemann, Antti Knowles 
and Anthony Metcalfe for fruitful discussions. During this work, AH was supported by the 
grant KAW 2010.0063 from the Knut and Alice Wallenberg Eoundation. The work of WH and 
JN was partially supported by the program “modeles numeriques” of the Erench Agence Na¬ 
tional de la Recherche under the grant ANR-12-MONU-0003 (project DIONISOS). Support 
of Labex BEZOUT from Universite Paris Est is also acknowledged. 

2 Statement of the main results 

2.1 The matrix model and assumptions 

The random matrix model of interest here is the N x N matrix 

( 2 . 1 ) 

where X^r is a A x n matrix with independent and identically distributed (i.i.d.) entries with 
zero mean and unit variance, and Sjv is a n x n deterministic positive definite Hermitian 
matrix. The random matrix Mjv thus has N non-negative eigenvalues, but which may be of 
different nature: The smallest N — min(n, N) eigenvalues are deterministic and all equal to 
zero, whereas the other min(n, N) eigenvalues are random. The problem is then to describe 
the asymptotic behaviour of the random eigenvalues of Mjv, as the size of the matrix grows 
to infinity. As for the asymptotic regime of interest, we let both the number of rows and 
columns of X^r grow to infinity at the same speed: We assume n = n{N) and n, A —)• oo so 
that 

Tl 

lim ^ = 7 e (0,oo) . (2.2) 

N^oo iV 
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This regime will be simply referred to as —)• oo in the sequel. 

Let us mention that the n x n random covariance matrix 


Miv = > 

which is also under consideration, has exactly the same random eigenvalues as Mat, and hence 
results on the random eigenvalues can be carried out from one model to the other immediately. 

Our first assumption is that the entries of are complex Gaussian. As we shall state 
later on, this assumption is fundamental for our local eigenvalue behaviour analysis, but not 
for our results on the limiting density behaviour, see Remark 2.3. 

Assumption 1. The entries of X^r are i.i.d. standard complex Gaussian random variables. 
Gonsidering now the matrix we denote by 0 < Ai ^ ^ An its eigenvalues and by 

1 

- (2.3) 

i=i 

its spectral measure. We also make the following assumption. 

Assumption 2. 

(a) For N large enough, the eigenvalues of Sjv stay in a compact subset of (0,+oo) inde¬ 
pendent of A, i.e. 

0 < liminfAi , sup An < -boo. (2-4) 

N^OO ]\J 

(b) The measure vn weakly converges towards a limiting probability measure v as N ^ oo, 
namely 

n 

(2-5) 

i=i 

for every bounded and continuous function /. 

Again, Assumption 2(a) is necessary for our results on the local eigenvalue behaviour, but 
our results on the limiting density behaviour require a weaker assumption, see Remark 2.3. 
We now turn to the description of the asymptotic eigenvalue distribution. 


2.2 Limiting eigenvalue distribution 

Consider the empirical distribution of the eigenvalues (xi) of Mjv, namely 

1 ^ 

i=l 

Since the seminal work of Marcenko and Pastur [25], it is known that this measure almost 
surely (a.s.) converges weakly towards a limiting probability measure fj, with compact support, 
provided that Assumption 2 holds true: 


1 ^ 


a.s. 

-^ 

N^oo 



( 2 . 6 ) 
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for every bounded and continuous function /. As a probability measure, /x is characterized 
by its Cauchy-Stieltjes transform, which is the holomorphic function defined by 

m{z) = f -- fJ,{dX), z G C+ = {z G C : Jm(z) > O} . (2-7) 

J z — X 

Marcenko and Pastur proved that m{z) is the unique solution m G C_ = {z G C : 3m{z) < 0} 
of the fixed-point equation 

where we recall that 7 has been introduced in (2.2) and v is the weak limit of see (2.5). 
Thanks to this equation, Silverstein and Choi then showed in [31] that /x({0}) = (1 — 7 )"*" and 
= m{x) exists for every x G M* = M — {0}. Consequently, the function m{z) 
can be continuously extended to C+ U M* and, furthermore, ^ has a density on (0, 00 ) given 
by 

p{x) =-Jm (m(x)). (2.9) 

vr 

We therefore have the representation 

/x(dx) = (1 — 7 )''' 5o + p{x)dx . ( 2 - 10 ) 

They also obtained that p{x) is real analytic wherever it is positive, and they moreover 
characterized the (compact) support Supp(/ 9 ) of the measure p{x)dx by building on ideas 
from [25]. More specifically, one can see that the function m{z) has an explicit inverse (for 
the composition law) on m(C+) given by 

g{m) = —+-/[ ^ ;/(dA) . ( 2 . 11 ) 

m 7 1 — mX 

If we introduce the open subset of the real line 

72 = {x G M : x / 0, x~^ ^ Supp(ix)} , (2-12) 

then the map g analytically extends to C+ U C_ U D. It was shown in [31] that 

M — Supp(p) = {s'(m) : m G D, g'{m) < O}. (2-13) 

Equipped with the definitions of m, g and D, we are now able to state our results con¬ 
cerning the behaviours of the limiting density p{x) near a cusp point or at the hard edge. 

2.3 Density behaviour near a cusp point 

As stated in the introduction, we define a cusp point 0 as an interior point where the density 
vanishes, namely 0 G int(Supp(/?)) such that /?(a) = 0. In particular, c = m(a) G M by virtue 
of (2.9). Our first result states that the density p{x) behaves like a cube root near a cusp 
point, provided that c G D. 
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Figure 1: Plot of p{x) with parameters 7 ~ 0.336 and v = 0.75i + 0 . 3(53 with a cusp point a. 


Proposition 1. Let a G int(Supp(p)) he such that p(a) = 0, and assume that c = m{a) G D. 
Then we have 

5(0 = 9 '(c) =/(O = > 0 • 

Moreover, 

= l^-ar^"(l + o(l)) , (2.14) 

In particular, there exists p > 0 such that for every x G (a — 7 , a + 7 ) \ {a}, we have p{x) > 0. 

Remark 2.1. In the forthcoming local analysis for the random eigenvalues near a cusp point, 
we shall focus on cusp points a’s satisfying a regularity condition. This extra assumption 
automatically yields that m(a) G D, see Remark 2.5. 

Conversely, we have the following result. 

Proposition 2. If c G D satisfies g'{L) = g''{c) = 0, then a = g{(.) belongs to int(Supp(p)) 
and p{a) = 0. In particular, g^^\c) > 0 and p{x) satisfies (2.14). 

We prove Propositions 1 and 2 in sections 3.1 and 3.2 respectively. Their proofs are based 
on the fact that there is a strong relation between the property that a = g'(c) is a cusp point 
and the local behaviour of g near c. For an illustration of these propositions, we refer to Figure 
2 where we displayed the graph of the map g associated with the density from Figure 1. 

We now turn to the hard edge setting. 
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Figure 2: Plot of function x i—?• g{x) on D for 7 ~ 0.336 and u = 0.7(5i + O. 363 . The vertical 
dotted lines are g^s asymptotes at x = 1/3 and x = 1. The thick segment on the vertical axis 
represents Supp(/?) who is delimited by the local extrema of g (see Eq. (2.13)). The point a 
is a cusp point. 

2.4 Density behaviour near the hard edge 

As usual in random matrix theory, the hard edge refers here to the origin when it belong 
to the limiting spectrum. In general, the limiting eigenvalue distribution may not display a 
hard edge, like in Figure 1. In fact, this is always the case when 7 7 ^ 1, see [17, Proposition 
2.4 (a),(c)]. Our next result states that there is a hard edge, namely 0 G Supp(/?), if and only 
if 7 = 1, and that in this case p{x) blows up like an inverse square root at the origin. We 
furthermore relate the presence of a hard edge to the behaviour of g near 00 . More precisely, 
since Supp(z/) C (0, 00 ) by Assumption 2, one can see from the definition (2.11) of g that the 
map g{l/z) is holomorphic at the origin. Thus, we have the analytic expansion as z —0 

gC^/z) = g{oo) + g'{oo)z + ^ ^ - (2.15) 

Clearly, g{oo) = 0 and the coefficients g'{oo) and g”{oo) are respectively given by the first 
and second derivative of the map z 1 -^ g{l/z) evaluated at z = 0 . 

Proposition 3. The following three assertions are equivalent: 

i) 0 G Supp(p) 

ii) 7 = 1 

iii) g'{oo) = 0 
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Moreover, we have g”{oo) = —27 f A ^z^(dA) < 0 and, if one of these assertions is satisfied, 
then 

1 / 2 A 

p{x) = -( —a;"^/ 2 (l + o(l)) , x-^ 0 +. (2.16) 

TT \-g”[oo)J 

Proposition 3 is proven in Section 3.3. 

Remark 2.2. There is an analogous statement for any left edge 0 > 0 of the spectrum 
satisfying c = m(a) G D which follows from [31]; see also [18, Section 2] for further information. 
Indeed, in this case we have 

5(0 = a, g'{c) = 0, g"{c) < 0, 

and furthermore, as x ^ a+, 

1 / \ 

^ (=7w) h'-“) + 

By analogy with this equation, the preimage c £ D corresponding to the hard edge is c = 00 . 
The fact that it actually belongs to D follows from Assumption 2(a). 

Remark 2.3. As we shall see in Section 3, proofs of Propositions 1, 2 and 3 only rely on 
the properties of the limiting eigenvalue distribution p, which do not depend on whether the 
entries of Xjv are Gaussian or not. More precisely, the exact assumptions required for these 
propositions are that the entries of X^r are i.i.d centered random variables with variance one. 
Assumption 2-(b), and that /.i({0}) = 0 (which follows from Assumption 2-(a)), see [25, 31]. 

2.5 The Pearcey kernel and fluctuations near a cusp point 

Our next result essentially states that the random eigenvalues of M^r, properly scaled near 
a regular cusp point, asymptotically behave like the determinantal point process associated 
with the Pearcey kernel, provided that an extra condition on a speed parameter is satisfied. 

In order to state this result, we first introduce this limiting point process. Next, we define 
what we mean by regular, and provide the existence of appropriate scaling parameters. After 
that, we finally state our result for the fluctuations near a cusp point. 


2.5.1 Determinantal point processes 


A point process on M (or in a subset therein), namely a probability distribution P over the 
locally finite discrete subsets (y*) of M, is determinantal if there exists an appropriate kernel 
K(x, y) : M X M —?• M which characterizes the correlation functions in the following way: For 
every k ^ 1 and every compactly supported Borel function $ : —)> M, we have 


E 


«>(yii,...,y4 


1 k 


- ■ ■ ,yk)det K{yi,yj) dyi • • • dy^ . 

L J i,j = l 


(2.17) 


In particular, the gap probabilities can be expressed as Fredholm determinants. Namely, 
given any interval J C M, the probability that the point process avoids J reads 


nj = 0 j =1 + ^ 


fc=i 


(- 1 )^ 

kl 


det 


'J JJ 


1 k 


K(yi,yj) dyi---dyfc, 
J i,j=i 


(2.18) 
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and the right hand side is the Fredholm determinant det(/ — K)^ 2 (j) of the integral operator 
acting on Lp‘{J) with kernel K(x, y), provided that it makes sense. For instance, if one assumes 
that J is a compact interval, which is enough for the purpose of this work, then det(I—K)£, 2 (j) 
is well-defined and finite as soon as ||K||y = sup^, |K(x, y)| < oo. Moreover, the map 
K(x,y) I—>■ det(/ — K)/^ 2 (j) is Lipschitz with respect to || • ||j when restricted to the kernels 
satisfying ||K|| j < oo (see e.g. [2, Lemma 3.4.5]). We refer the reader to [19, 21, 2] for further 
information on determinantal point processes. 


2.5.2 The Pearcey kernel 

Given any r G M, consider the Pearcey-like integral functions 




1 

2z7r 



^XZ-T^ 


4 dz, 


iiv) = ^ 


f 


^-yw+T ^— 


4 dw, 


where the contour S has two non-intersecting components, one which goes from 
to whereas the other one goes from to More precisely, we 


parametrize here this contour by 

S = |±e*® : 0 G [—7r/4,7r/4]| U t G (—oo, —1] U [1,-|-oo)| , (2-19) 

with the orientation as shown in Figure 3. 




Figure 3: The contour S. 

It follows from their definitions that the functions (p and ip satisfy the respective differential 
equations 

- Tp'ix) + X(j){x) = 0, Ip”'{y) - Tip'{y) - yip{y) = 0. 

The Pearcey kernel is then defined for x, y G M by 

(p"{x)ip{y) - (j)'{x)ip'{y) -h (j){x)ip"{y) - T(p{x)ip{y) 




x-y 


( 2 . 20 ) 


see for instance [8] . One can alternatively represent this kernel as a double contour integral 


Kpj {x, y) = 


{2m) 


dz 


dw 


gXZ 2 


2 4 2 4 

-- ■ - 


w — z 


( 2 . 21 ) 
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from which one can easily see the symmetry = K.p^{—x,—y) by performing the 

changes of variables z —z and w ^ —w. 

The Pearcey kernel first appeared in the works of Brezin and Hikami [7, 8 ] when studying 
the eigenvalues of a specific additive perturbation of a GUE random matrix near a cusp point. 
Subsequent generalizations have been considered by Tracy and Widom [34], and a Riemann- 
Hilbert analysis has been performed by Bleher and Kuijlaars [5] as well. This kernel also 
arises in the description of random combinatorial models, such as certain plane partitions 
[27]. Furthermore, it has been established that the gap probabilities for the associated point 
process satisfy differential equations. For instance, logdet(/ —satisfies PDFs with 
respect to the variables s, t and r, see [34, 4, 1], which should be compared to the connection 
between the Tracy-Widom distribution and the Painleve II equation. 


2.5.3 The regularity condition 

We start with the following definition. 

Definition 2.4. A cusp point o is regular if c = m(a) satisfies 


n 

lim inf min 

N^OO j = l 



( 2 . 22 ) 


The regularity condition (2.22) has been considered in [17] when dealing with soft edges, to 
ensure the appearance of the Tracy-Widom distribution. Similarly here, as we soon shall see, 
this condition enables the Pearcey kernel to arise at a cusp point. Moreover, the behaviour of 
p{x) at such cusp points is well described by Proposition 1, as explained in the next remark. 


Remark 2.5. If a is a regular cusp point, then it follows from the weak convergence ^ v 
and the definition of D that necessarily c = m(a) G D. Thus, o satisfies the hypothesis of 
Proposition 1. In particular, we have a = g{c) and p{x) behaves like a cube root near a. 

Finally, the regularity assumption yields the existence of natural scaling parameters for 
the eigenvalue local asymptotics. Consider the counterpart of the map g introduced in (2.11) 
after replacing zz by see (2.3), and 7 by n/N, namely 

9n{z) = ^ ^ y ^NidX) . (2.23) 


The map g'Ar is the inverse Cauchy transform of a probability measure pN usually referred to 
as deterministic equivalent for the random eigenvalues distribution of Mv, see [18, Section 
3.2]. According to Section 2.2, pj\f has a decomposition of the form (2.10). In particular, it 
has a density pN on ( 0 , + 00 ) which is analytic wherever it is positive. 

The next proposition provides an appropriate sequence of finite-A^ approximations of c, 
which we will use in the definition of the scaling parameters. 


Proposition 4. Let a = g{c) be a regular cusp point. Then there exists a sequence (cjv) of 
real numbers, unique up to a finite number of terms, converging to c and such that for every 
N large enough, we have ci\f € D and 


lim gNicN) = gic) , lim g'^icN) = 0 , 
N^oo N^oo 


dNi^N) — 0 


and 


> 0 . 
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This proposition is the counterpart of [17, Proposition 2.7], with a similar proof. Let us 
only provide a sketch here: Combined with Montel’s theorem, the regularity condition ensures 
that converge uniformly to on a neighbourhood of c for every /c ^ 0. The proposition 
then follows by applying Hurwitz’s theorem to g'^ since c is a simple root for g", according 
to Proposition 2.14. 

Let us emphasize that there is however an important difference with regular soft edges as 
described in [17], where it was shown that if a = g{c) is a regular soft edge for g, then 5 'Ar(cAr) 
is a soft edge for gj\f. 

Remark 2.6. A regular cusp point may not be the limit of finite-A^ cusp points. More 
precisely, if o = g{c) is a regular cusp point, then in particular g\c) = g”{c) = 0. However, 
this only ensures the existence of a sequence (cat) such that g'^{cN) = 0. A priori, g']y(cAr) 

converges to zero as Ai —)• oo but might not be equal to zero. In fact, it is not hard to show 

we have the following alternatives: 

• if g'j\f{cN) = 0, then gwi^^N) is a cusp point for /ijv; 

• if 5 (]Y(cAr) > 0, then the density pat is positive in a vicinity of PAf(cAf); 

• if p]y(cAr) < 0, then gNi<^N) does not belong to the support of g^- 

We are finally in position to state our result concerning the eigenvalue behaviour at a 
regular cusp point. 


2.5.4 Fluctuations around a cusp point 

Thanks to Assumption 1, the random eigenvalues (xi) of Mtv form a determinantal point 
process with respect to a kernel K]sf{x,y), see [3, 28]. An explicit formula for this kernel is 
provided in Section 4. The main result of this section is the local uniform convergence of this 
kernel, properly scaled, towards the Pearcey kernel. 


Theorem 5. Let a = p(c) be a regular cusp point. Let (cat) be the sequence associated to c 
coming from Proposition 4- Assume moreover that the following decay assumption holds true: 
There exists k G M such that 


^5'jv(CAf) —-^ ■ 

N^oo 


(2.24) 


Set 


1/4 


aiv = fl'Af(CAf) , 


cjat = 




T = —K 


g(3)(c 


1/2 


(2.25) 


so that Oat —)• a and —>• (6/p(^)(c))^/'^ > 0 as N ^ oo by Proposition 4- Then, we have 


uniformly for x,y in compact subsets o/M. 

This result was obtained by Mo [26] in the special case where the matrix 'Sj\[ has exactly 
two distinct eigenvalues (each with multiplicities proportional Ai), by means of a Riemann- 
Hilbert asymptotic analysis. 
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Notice that if (xj) is determinantal with kernel K 7 v(x, 7 /), then — ajv)) is 

determinantal with the kernel given by the left hand side of (2.26). Thus, having in mind 
Section 2.5.1, a direct consequence of this theorem is the convergence of the compact gap 
probabilities. 

Corollary 2.7. Under the assumptions of Theorem 5, we have for every G M, 

P ((iV3/V,v(xi - UTv)) n [s,t] = 0 ) det(I - • (2-27) 

We now make a few comments on the assumption (2.24). 

Remark 2.8. The decay assumption (2.24) roughly states that the cusp point a = g{c) 
appears fast enough. More precisely, for a cusp point a = g{c), one has in particular g'{c) = 0. 

• If K > 0, then g'jsfi^N) > 0 for large N. According to Remark 2.6, the density pjy is 
positive near gwi^^N) and will converge to zero to asymptotically give birth to a cusp 
point. The family of densities pN display a sharp non-negative minimum at S'Ar(cAr) 
converging to zero which may be thought of as the erosion of a valley, see the thin curve 
in Figure 4. 

• If K < 0, then g']\i{cN) < 0 for large N and the density pN vanishes in a vicinity of 
fl'Af(cAf)- However, this interval will shrink and asymptotically disappear. Thus, two 
connected components of the support of pjsf move towards one another (moving cliffs), 
see the dotted curve in Figure 4. 

The assumption (2.24) is an indication on the speed at which the bottom of the valley reaches 
zero (k > 0), or at which the two cliffs approach one another (k < 0). See [18] for a more 
in-depth discussion. 


Remark 2.9. When the assumption (2.24) is not satisfied, namely when g'j^{c]sf) goes to zero 
as N 00 (which is always true by Proposition 4) slow enough so that \/lVg'j^{cN) diverges 
to plus or minus infinity, we do not expect the Pearcey kernel to arise. See [18] for further 
discussion. 

The proof of Theorem 5 is provided in Section 4. 

2.6 Asymptotic expansion at the hard edge 

Our last result concerns the behaviour of the smallest random eigenvalue of when the 
hard edge is present. Recall that n/N —)• 7 . By Proposition 3, the limiting density displays 
a hard edge if and only if 7 = 1. With this respect, we restrict ourselves here to the case 
n = N + a where a G Z is fixed and does not depend on N. 


2.6.1 The Bessel kernel 

The Bessel function of the first kind Ja with parameter a G Z is defined by 


j _til!_ 

Z^^!r(n + a + l) I2J ’ 

n=0 ^ ' 


X > 0 , 


(2.28) 
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Figure 4: Zoom of the density of /tat near the cusp point a. The thick curve is the density of 
^ in the framework of Figure 1. The thin curve (resp. the dotted curve) is the density of 
when > 0 (resp. VNg'^{cN) < 0). 


with the convention that, when a < 0, the first |q;| terms in the series vanish (since the 
Gamma function F has simple poles on the non-positive integers). 

The Bessel kernel is then defined for x,y > 0 by 




y/yJa{Vx)j'^{y/y) - y/xJ'„{y/x)Ja{^) 
2{x - y) 


(2.29) 


One can alternatively express it as a double contour integral, 


K£^(a:,2/) = 


1 


(2i7r)^ 


|2|=r ^ 


drc 1 


— e 


(2.30) 


W z-w 

where 0 < r < R and the contours are simple and oriented counterclockwise, see for instance 
[17, Lemma 6.2]. We set for convenience 


Fa{s) = det (/ - k£^)^ 2 ( 0 ,^), s > 0, (2.31) 

where the right hand side stands for the Fredholm determinant of the restriction to L^(0,s) 
of the integral operator According to (2.18), Fa{s) is the probability that the smallest 

particle of the determinantal point process associated with the Bessel kernel is larger than 
s. Tracy and Widom [33] established that certain simple transformations of logF’„(s) satisfy 
Painleve equations (Painleve III and Painleve V are involved). 


2.6.2 Correction for the smallest eigenvalue’s fluctuations 


We denote by Xmin the smallest random eigenvalue of Mv, namely 


^min 


xi if a ^ 0, 

xi-a if a < 0. 


(2.32) 
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Our last result is stated as follows. 


Theorem 6. Assume n = N + a where a € Z is fixed and does not depend on N. Set 


4^1 8^1 

■ ' i=l * 


j=l ^ 


so that 


O'N 


-)> 4 / —^— > 0 and Cn 


N^oo J X " N^oo 

by Assumption 2. Then, for every s > 0, we have as N ^ oo, 


z^(dA) 


> 0 


(2.33) 


P^^iY^fJATXmin ^ ^ ■ (2.34) 


The convergence towards Fa(s) has been first observed by Forrester [15] when is the 
identity. As for the general case, it has been established by the authors in [17]. An explicit 
formula for the 1/A^-correction term was conjectured when S^v is the identity by Edelman, 
Guionnet and Peche [12], a conjecture proved true soon after by Perret and Schehr [29] and 
Bornemann [6], with different techniques. We thus generalize this formula to the general S^r 
case. The strategy of the proof is rather similar to Bornemann’s one: It relies on an identity 
involving the resolvent of obtained by Tracy and Widom, although we cannot rely on 
existing estimates for the kernel in this general setting. 


Remark 2.10. In fact, as we shall see in the proof of Theorem 6 (see Remark 5.2), our 
method easily yields for every s > 0 an expansion of the form 


P^A^ (Tat Xmin ^ — Fa{s) + + 0 f (2.35) 

£=1 ^ ^ ^ 

for every L ^ 1 as A —)• oo. Although we are able to provide a close formula for the 
coefficient (17^1(5) (as stated in Theorem 6) thanks to a formula due to Tracy and Widom, 
to the best of our knowledge the next order coefficients do not seem to benefit from such a 
simple representation. 


We prove Theorem 6 in Section 5. 


3 Proofs of the limiting density behaviours 

This section is devoted to the proofs of Proposition 1, 2 and 3. 

We first recall a few facts stated in Section 2.2 that we shall use in the forthcoming 
proofs: The map m is the Cauchy-Stieltjes transform (2.7); it is analytic on C+ and extends 
continuously to C+UM*. Moreover, iJm(m(x)) = —7rp(x) for every x G M*. The map g defined 
in (2.11) is analytic on C+ U C_ U D. In particular, g' has isolated zeroes on D. Moreover, 
after noticing that m(C_|_) C C_, we have the identity g{m{z)) = z for every z G C+. 

We start with a simple but useful fact, which follows by taking the limit z G C+ —)• x in 
the previous identity and using the continuity of m and g on their respective domains: 

Lemma 3.1. If x (z C+ U M* is such that m{x) G C_ U D, then g{m[x)) = x. 
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We will also use the following property. 

Lemma 3.2. //a G Supp(/o) satisfies p{a) = 0 and c = m{a) G D, then g'{c) = 0. 

Proof. Consider the map 

G'(m) =+7 / ——^z^(dA) , m G C_ U D, (3.1) 

Imp J |1 — Amp 

and notice it is continuous on its definition domain by dominated convergence. It follows 
from the dehnition (2.11) of g that G{m) = g'{m) when m D, and moreover that for every 
m G C_ we have the identity 


iJm(g'(m)) = 3m{m)G{m). 

By using the fact that g{m{z)) = z on C+, we thus obtain for every 2 G C+, 

Jm(z) = '3m{m{z))G{m{z)). (3.2) 

If X G M* is such that p{x) > 0, then by letting z G C+ —)• x in (3.2) we see that necessarily 
G{m{x)) = 0 because lim 2 _>.a; Jm(m(z)) = Jm(m(x)) = —7rp{x) < 0. Now, since 0 G Supp(p) 
by assumption, there exists a sequence (x^) such that x^ —)• a as /c —)■ 00 and p{xk) > 0, and 
hence G{m{xk)) = 0. Since m{xk) —>■ m(a) = c and c £ D, this yields gfc) = G(c) = 0. □ 

3.1 Density behaviour near a cusp point: Proof of Proposition 1 

We now turn to the proof of the first proposition. 

Proof. Assume that 0 G int(Supp(/?)) and /9(a) = 0. Set c = m(a) and assume moreover that 
c £ D. Thus, the facts that g{c) = a and g'{c) = 0 directly follows from Lemma 3.1 and 
Lemma 3.2. 

First, we prove that g"{c) = 0. To do so, we show that < 7 ' > 0 on (c — 7 , c + 7 ) — {c} for 
some 7 > 0. Since g'{c) = 0 this would indeed yield that c is a local extremum for g'. We 
proceed by contradiction: Assume there exists a sequence (xk) in — {c} such that x^ —>■ c 
and g'{xk) ^ 0. Since g' has isolated zeroes on D, necessarily g'{xfi) < 0 for every k large 
enough. It then follows from (2.13) that g{xk) G M — Supp(/ 9 ) and, since g{xk) —)■ fl'(c) = a, 
this contradicts the assumption that 0 G int(Supp(p)). 

Next, we similarly show that there exists 7 > 0 such that 

X G (a — 7 , a + 7 ) — {a} ^ p{x) > t) . (3.3) 

We will use (3.3) later on in this proof. Assume there exists a sequence (x^) in M* — {a} 
such that Xfc —>• 0 and p{xk) = 0. Since a G int(Supp(/ 9 )) and m(a) G H by assumption, then 
Xk £ Supp(p) and m{xk) £ D for every k large enough. Moreover, we have m{xk) —>■ fn{a), 
but since Lemma 3.2 then yields g'{m{xk)) = 0, this contradicts that g' has isolated zeroes 
on D and (3.3) follows. 

Now, we show that g^^\c) > 0 by direct computation. Recalling the definition (2.11) of 
g, the equation g”{c) = 0 reads 


15 





As a consequence, we obtain 





A4 

(1 - Ac)4 


z^(dA) 


We finally turn to the proof of the cube root behaviour (2.14). Since g'{c) = g"{c) = 0 
and g^^\c) > 0, there exists an analytic map tp defined on a complex neighbourhood of c 
such that g{m) — a = p{m)^ and p>'{c) / 0, see e.g. [30, Theorem 10.32]. In particular, 

we have p'{c)^ = g^^\c)/Q = C > 0. Moreover, the inverse function theorem yields that 
p has a local inverse defined on a neighbourhood of zero, such that (/?“^(0) = c and 
(0) = l/ip'{c). If \x — a| is small enough, then p{x) > 0 by (3.3), hence m{x) G C_, 
and we have g{m(x)) = x by Lemma 3.1. In particular, we have 


(p{m{x))^ = g{m{x)) — a = x — a 


and, by taking the cube root (principal determination), applying p and performing a 
Taylor expansion, we obtain 

x-ay/3 

) 


+ o{x — a)^/^ 


(3.4) 



where j is an undetermined cube root of unity. Since m{x) G C_, necessarily j = exp(2f7r/3) 
if a: < a and j = exp(—2i7r/3) if x > a. Finally, (2.14) follows by taking the imaginary part 
in (3.4). 

Proof of Proposition 1 is therefore complete. 

□ 


3.2 Identification of a cusp point: Proof of Proposition 2 

The main part of the proof consists in showing the following lemma. 

Lemma 3.3. Let c € D such that g'{c) = 0. Then g(c) ^ 0 and m{g{c)) = c. 

Equipped with Lemma 3.3, let us first show how the proposition follows. 

Proof of Proposition 2. Assume that c £ D satisfies g'{c) = g”{c) = 0 and set a = ^(c). 
In particular, a G Supp(p) by (2.13), and we just have to show that p(a) = 0 and a ^ 
dSupp(/9). We know from Lemma 3.3 that a 7 ^ 0 and m(a) = c. As a consequence, p(a) = 
—Jm(m(a))/ 7 r = 0. Finally, since a = g{c) with c £ D and p(o) = 0, then [31, Theorem 5.2] 
shows the condition a G 9Supp(/9) requires g”{c) / 0, which is not possible by assumption. 

□ 


We now prove the lemma. 
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Proof of Lemma 3.3. For every m € D we have 


g{m) + mg{m) = j J ^(dA) / 0 , 

and because g'{c) = 0 by assumption, by taking m = c we see that necessarily g{c) 0. Let 
us now prove that m{g{c)) = c. Introduce for convenience the map 


$(z) 


1-z ■ 


By combining the fixed point equation (2.8) for m{z) with z G C+ and that 


1 

c 


g{c)-g{c) + 


1 

c 


9 ( 0-7 j 


A 

1 — cA 


z/(dA) , 


we obtain 


m(z) — c = cm(z) (- 

' c m{z) 


= cm{z){g{c) — z) +'y{m{z) — c) / <h(Am(z))<I>(Ac) z/(dA) 


By reorganizing this equation as 


{m{z) — c) ^1 — 7 J <h(Am(z)) $(Ac) z^(dA)^ = {g{c) — z)m{z)c , (3.5) 


we see the lemma would follow by taking the limit z G C-|- —)• g{c) G M* assuming that 


1 — 7 y <k(Am(z)) <I>(Ac) i^(dA) 


^ C'|m(z) — cp , z G C+ —)• fl'(c), 


(3.6) 


for some constant C > 0. To show (3.6), we start by writing 


1 —'J J <I>(Am(z))‘k(Ac) i/(dA) 

^ 1 —'J J |<k(Am(z)) <l>(Ac)| i/(dA) (3.7) 

|<I>(Am(z))|^ i/(dA) + J |<h(Ac)|^ i/(dA) — y |$(Am(z)) — <I>(Ac)|^ z^(dA) 


= 1 -^ 


Recalling the definition (3.1) of G, we see from (3.2) that G{m(z)) < 0 on C+, and this yields 

limsup 7 j |<I>(Am(z))|^ zz(dA) ^ 1 . (3.8) 

Similarly, the equation G(c) = g^c) = 0 reads 

7 y |$(Ac)|^z/(dA) = 1 . (3.9) 
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Finally, we have 


J |^>(Am( 2 ;)) — <h(Ac)|^i/(dA) = \m{z) — c|^ 

^ \m{z) — c|^ 


A 


(1 — Xm{z)){l — Ac) 
A 


(1 +A|m(z)|)(l + A|c|) 
and moreover, since m is continnous on C+ U M* and g{c) G M*, 

A 2 


i^(dA) 

u{dX) (3.10) 


lim inf 
2eC+-5>g(c) , 


i^(dA) > 0 . 


(3.11) 


(1 + A|m( 2 ;)|)(l + A|c|) 

The lower bonnd (3.6) then follows by combining (3.7)-(3.11) and the proof is complete. □ 


3.3 Density behaviour near the hard edge: Proof of Proposition 3 

The key to study the hard edge is to study the map g near infinity, which is holomorphic 
there. More precisely, we have the analytic expansion as z —)■ 0, 


g{l/z) = z-'yz j 


\ — zjX 


z^(dA) 


= (1 - 7)z - (7 y A-^z^(dA)^ ^^ + • • • (3-12) 


The reason why one should do so is that \m{x)\ goes to 00 as x \ 0: 

Lemma 3.4. Assume 0 G Supp(/?). Then, there exists 7 > 0 sueh that p{x) > 0 on (0,r/). 
Moreover, as x > 0 decreases to zero, we have |m(x)| —)> + 00 . 


Equipped with this lemma, we first provide a proof of the proposition. 

Proof of Proposition 3. A comparison between (2.15) and (3.12) readily yields g'{oo) = 1 — 7 
and ^''(oo) = —27 f X~^v(dX) < 0. In particular, the equivalence between ii) and Hi) is 
obvious. That i) =i> ii) follows from [17, Proposition 2.4(a) (c)]. We show ii) ^ i) by 
contradiction: Assume 7 = 1 and that there exists 0 < oq < minSupp(/ 3 ). Since 7 = 1, by 
(2.10) we have g{dx) = p{x)dx and thus cq = m{ao) = f (oq — x)~^p(x)dx < 0 . A close look 
at the definition (2.11) of g shows that <7 < 0 on (—oo,0) (see also [17, Proposition 2.4(b)]), 
and thus g(co) < 0. On the other hand, cq € D since cq < 0 and hence Lemma 3.1 then yields 
0 < ao = g(co), which is a contradiction. 

We now prove the inverse square root behaviour (2.16). Since g'(oo) = 0 and g"(oo) < 0, 
there exists an analytic map (p defined on a complex neighbourhood of the origin such that 
g(l/z) = ip(z)^ and <^^(0) 7 ^ 0, see [30, Theorem 10.32]. In particular, we have ip'(O)^ = 

g"{oo)/2 = —C < 0. The map (p is uniquely defined up to a sign, and thus p>'{d) = 
Moreover, ip has a local inverse ip~^, defined on a neighbourhood of zero, so that (/j“^( 0) = 0 
and (<F~^) (0) = \/ip'{d). For every x > 0 small enough. Lemma 3.4 yields that m(x) G C_, 
and thus g{m{x)) = x by Lemma 3.1, and moreover that l/m(x) lies in the definition domain 
of (/J. As a consequence, 

ip{l/m{x))‘^ = g{m{x)) = x 
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and, by taking the square root (principal determination), applying ip and performing a 
Taylor expansion, we obtain 

Since m[x) G C_, the undetermined sign has to be a minus sign, and (2.16) follows by taking 
the inverse and then the imaginary part in (3.13). 

□ 


We finally turn to the proof of the lemma. 

Proof of Lemma 3. 4- Let z = x + iy with x,y > 0 and assume that p{x) > 0. The fixed point 
equation (2.8) yields 

= -91e(m(z))G(m(z)) + 7 y zy(dA), 

where G is defined as in (3.1). When y —>■ 0, since p{x) > 0, we have m{z) —)> m{x) G C_ and 
moreover G{m{z)) —)• 0 for the same reason as in proof of Lemma 3.2. Thus, we obtain. 

We now show that p{x) > 0 for every x G (0, p) for some y > 0 by contradiction: Assume 
there exists a sequence (x^) such that x^ > 0, x^ —)> 0 as A: —>■ oo, and p{xk) = 0. Since 
0 G Supp(/9), without loss of generality one can assume that Xk G Supp(/3). As a consequence, 
using moreover that m is continuous on M*, one can construct a sequence {yk) satisfying 
yfc > 0, yfc —>• 0 as A: —)■ oo, p{xk) > 0, and |m(xfc) — m(yfc)| ^ 1- Because Assumption 
2 yields minSupp(zA) > 0, we see from (3.14) that necessarily \m{yk)\ —>■ +oo, and hence 
|m(xfc)| —>■ + 00 , as A: —)> oo. In particular, m{xk) G H for A; large enough and Lemma 3.2 
then yields g'{m{xk)) = 0. Thus, (l/m(xA;)) is a vanishing sequence as A: —>• oo of zeroes of 
the derivative of g{l/z). But since g{l/z) is holomorphic near zero, and hence so does its 
derivative, this contradicts the isolated zero principle and our claim follows. 

Finally, since the identity (3.14) thus holds true on (0,r/), that |m(x)| —)> +oo follows by 
letting X > 0 decrease to zero in this identity. □ 


4 Fluctuations around a cusp point: Proof of Theorem 5 


In this Section, we study the local behaviour of the random eigenvalues (xj) of M^r near a 
regular cusp point and prove Theorem 5. Our asymptotic analysis is based on that, thanks 
to Assumption 1, the random eigenvalues form a determinantal point process with an explicit 
kernel K 7 v(x,y), see [3, 28]. The kernel has the following double contour integral formula 


K7v(x,y) 


N 

(2i7r)2 



Q-Nx{z-q)+Ny{w-q) 



(4.1) 


where the (7 G M is a free parameter and we recall that the Aj’s are the eigenvalues of 51 tv. T 
and 0 are disjoint closed contours such that T encloses all the A“^’s whereas 0 encloses the 
origin, see for instance Figure 5. 
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Convention: All the contours we consider are simple and oriented counterclockwise. 



Figure 5: The contours F and 0 in the definition of K 7 v(x,y) 


Remark 4.1. In [3, 28], it is assumed that q satisfies an extra restriction so that the associated 
operator Kjv is trace class on the semi-infinite intervals (s,-|-oo). Since we are here only 
interested in establishing a local uniform convergence for 'Kjq{x,y), such restrictions are not 
necessary. See also [17, Remark 4.3]. 

In the remaining of this section, we assume that the framework (assumptions and nota¬ 
tions) of Theorem 5 holds true. We then set 


(“" + JV^' + JV^ 


(4.2) 


and, in order to establish Theorem 5, focus on the proof of the uniform convergence 


sup 

x,yE[—s,s] 


^N{x,y) - K^^l{x,y) 


N^OO 


-> 0 


(4.3) 


for every fixed s > 0 . 


Notations and conventions: 


• We denote by B{z, p) the open disc in C with center z G C and radius p > 0. 

• By convention, we shall use at several instances a constant C which depends neither on 
N nor on x, y G [—s, s], but may depend on s, and whose exact value may change from 
one ligne to another. 


If a contour F is parametrized by 7 : / —)• F for some interval / C M, then for every map 
/i : F —?• C we set 

J h{z)\dz\ = j /i o 7 (t) | 7 '(t)|dt 

when it does make sense. In particular, /pjdzj is the length of a closed contour F. 


4.1 Preparation 

It follows from (4.1) and (4.2), by setting q = cat, that 


^N{x,y) 


ArV4 

(2i7r)2(TAr 





vi/4 

(TN 


y(w-CN) fN{z)-N f^iw) 


(4.4) 
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where we introduced the map 


fN{z) 


-aN{z - Cat) +log( 2 ;) 


1 

-^log(l-A,z) . 
i=i 


(4.5) 


Notice the functions exp(/Ar), iHe/Ar, and the derivatives (/c ^ 1) are well-defined on 
C \ {0, , A“^}. However, one needs to specify an appropriate determination of the 

complex logarithm in order to define fjsr properly. By Assumption 2 and the regularity 
condition (2.22), there exists e > 0 such that Aj^ G (0, -|-oo) \ B{c, e) for every 1 ^ j ^ n and 
every N large enough. As a consequence, if we introduce the compact set 


JC 


'■ f ^ 1 ' 

inf ^,sup — 
N An Af ^1. 


\H(c,e) U{0} 


(4.6) 


then on every simply connected open subset of C \ /C one can find a determination of the 
logarithm so that /at is well-defined and holomorphic there. 

Recalling the definition (2.23) of qm, an essential observation is that 


f'N{z) = gN{z) - aw ■ 


(4.7) 


As a consequence, we have for every k ^ 1, 

= g^^\z) . (4.8) 

In particular, since otv = fl'Af(cAf) by definition, the decay assumption (2.24) and Proposition 4 
provide 


o 

II 

(4.9) 

y^fNi'^N) —- ^ T , 

N^oo 

(4.10) 

= 0 , 

(4.11) 

> 0 . 

A'^oo 

(4.12) 


By performing a Taylor expansion of /at near cat in (4.4), one can already guess from 
(4.9)-(4.12) and a change of variables that the Pearcey kernel should appear in the large N 
limit, at least if one restricts the contours P and 0 to a neighbourhood of ca^ (see also [18] for 
a more detailed heuristic, which may serve as a guideline for the forthcoming proof). In a first 
step, we provide precise estimates in order to prove that claim, see Section 4.2. In a second 
step, we prove that the contribution coming from the pieces of contour away from cn are 
exponentially negligible, see Section 4.3. To do so, we establish the existence of appropriate 
contours in the same fashion as in [17]. This will enable us to conclude. 

Notice that one of the key arguments to assert the existence of appropriate contours is 
the maximum principle for subharmonic functions. This argument only appears in the proof 
of Lemma 4.9, which is similar to the proof of [17, Lemma 4.11] and hence omitted. The 
interested reader may refer to [17] for more details. 
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4.2 Step 1: Local analysis around cn 

We start with a quantitative Taylor expansion for /tv near ctv- 

Lemma 4.2. There exists po > 0 and A = A(po) > 0 independent of N such that for every N 
large enough, B{cn,Po) T C\JC and, whatever the analytic representation of f^ on B{cn, Po), 
we have for every z G B{c]\f,po) 

\fNiz) - fNiCN) - ^5'7v(Civ)(^ - (cat) - Cn)^\ ^ A|z - Ctv|^ 

In particular, since 5 ^(ctv) is real, for every z G B{cn,Po), 

|5te/Tv(^) - 9^e/Ar(cAr) - ^g'^{cN)^e{z - cat)^ - ^gP{cN)^^{z - ctv)^| ^ A|z - cat]'’. 

Proof Choose po such that B{c,2pq) C C \ /C. The convergence ctv —c then yields 
-S(cjv,/9o) C B{c,2pq) C C \ /C for every N large enough. By using (4.8)-(4.11) and per¬ 
forming a Taylor expansion for /tv around ctv, we obtain 


\fN{z) - /Af(CTv) - ^9'n{'^n){z - Ctv)^ “ (Civ) “ Ctv)^| 

= \fNiz) - /jv(ctv) - ^/^(CAr)(^ “ Ctv)^ “ ^fN\<^N){z - Ctv)"^ 

^ ^1 1 5 I .f (^) I 


provided that z G B{c]y,po) and N is sufficiently large. Moreover, since /j^^ = converges 
uniformly on B{c, 2po) to < 7 ^^^ which is bounded there, the existence of A = A(po) independent 
of N follows. □ 


From now, we let p > 0 be small enough so that 0 < p < po and 


We introduce the contours 


gff(cjv) _ ^ 1 

4! ^^24! 


(4.13) 


T+ = {cTv + iV-^/V': gG [-7r/4,7r/4]} , (4.14) 

Tl'x = jcAf — : t G [—p, — U jcAT + : tG[A“^/^,p]| , (4:.15) 

T; = {cTv-iV■^/V^ 0G [-7r/4,7r/4]} , (4.16) 

T- = |cTv + te“/^ : t G [-p,-A-V4]| y |ctv + : t G p]'j . (4.17) 


The orientations of these contours are specified by letting the parameters 9 and t increase in 
their definition domains. We also introduce 




Similarly, the orientation of contour is specified by letting t increase in its definition 
domain. 

We now establish the following estimate which essentially allows to replace /at by its 
Taylor expansion around cat in the double integral over the contours T* and H*. 


Lemma 4.3. For every s > 0, the following quantity 


ArV4 f , f , 1 -Ari/43. (^-^iv) I Afi/4^ (“-^iv) 

7 — / az aw -e <^n 

{2i'kY(Tn Jr, Js, w - z 

( / (3) / (3) ' 


(4.20) 


converges to zero as N ^ oo, uniformly in x,y € [—s,s]. 
Proof. Let s > 0 be fixed, and set for convenience 

Dn{z,w) = 




SO that it amounts to prove that 


sup 

3;,ye[-s,s] 


f 1 _Ari/4a.(£^+Ari/4j/Li^ ^ , , 

az / drc-e <^n D]\f{z,w) 

w - z 


vanishes as N ^ oo. 

First, by definition of the contours (4.14)-(4.19), we have the bound 


^ V2ivV4 


z — w 


for every z € T, and w G E,. As a consequence 
iVV4 


f . f . 1 -ArV4a; (^ '^Jv) _)_Afi/4^ (™ ^n) 

f dz / dw -e '^N Disf(z,w) 

r, Je, w- z 


^ \/2iV / \dz\ / |dtc| e 

«/ X* j .IT;* 


Next, we use the following elementary inequality, 

IgW _ _ g9te(r))|g(n-i;) _ 




Dn{z,w] 


C e 


tyie{v) 


E 

fe^i 


u — u 


k\ 


^ kt — u e 


— ,,|p5^c(t;)+|u-'u| 


(4.21) 


(4.22) 


(4.23) 


which holds true for every u,v G C. By specializing it to 

u = Nf^iz) - NfM{w) = N{fN{z) - fN{CN)) - N{fN{w) - /Ar(CAr)) , 

jQ'n^^n) , x2 I , ,^4 


V = N- 


2 V 41 V 

- c„ f - - w)‘ , 
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where z G T* and w G H*, so that Lemma 4.2 yields |tt — r;| ^ iVAdz — CAr|® + |t(; — cn\^) , we 
obtain 


Dn{z,w) 


^ AA(|z — ctvI^ + lii; — catI^) 

/ (^) 

^ ^N^^i^dKc(z-CN)^+N^!^-p^m,iz-CN)*+NA\z-CN\^ 

/ (3) 

^ We('»)-Civ)^-jV ^^4,‘^’ We('»)-Civ)"‘+jVA|«)-Civ|^ 


As a consequence, it follows 


A ^ / |dz| / Idle] e 

«/ X* t/ .n* 




Dn{z,w] 
(3), 


^ / AIz-cat 




L 


/ (^) 


X / e 


cr^ 


|dt(;| 


I 


/ (3) 

“1“ / C Cl2^ 

'T* 


/ (2) 
^i/4y iHU^-cjv) _jvgk^9^e(«)-c^)3-jV ^^4^;^^ 9re(io-c^)4+jVAh-cjvP 


X / A|t(; — CAr|^e 


|dt(;|. 

(4.24) 


We now handle the integrals over each piece of contour separately. 

First, consider the integrals over the contour H*. By definition tc G H* if and only if there 
exists t G [—p, p] such that tc = ca^ + it. Thus, 

yit{w — Cat) = 0, iRe{w — cat)^ = — y\.t{w — cn)^ = t*, \w — catI^ = |t|^- 

For N large enough this yields, together with (4.13), 

«/ a:.* 

/ (^) 


i-p 


,9']y(tN)j_2 ATii( 


€ 




p P-NP 

/ " 

I-p 

^^gAr('v)^2 T^Ti-j-4, 

/ dt 

I-p 


'-pA 


dt 


1 




Ai/4 y_^^i/4 

1 ,/jv4d!iv) 


AV4 /_ 


J —OO 


c 

AV4 


(4.25) 
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where we set for convenience 


1 g^^\c) 


For the last estimate, we have used the decay assumption \/]V g'j^{cN) —)• /t G M. Similarly, 


J N\w — cisfl^e 


/ (2) 


|drc| 


J — p 


'-P 

CP 


^ N / |tre^' 2 
J-p 


F r'" i 

VN J-pN^/i 


c 

Vn ■ 


(4.26) 


Next, we turn to the integrals over the contours and T“. By definition z G T((, 

resp. z G T“, if and only if there exists t G [N~^/^,p\ such that z = tM + resp. 

z = cn — so that in both cases we have 

|iHe(2; — Cat)! ^ t, lHe(2; — cat)^ = 0, 94e(2: — ctv)^ = —|^ —CAr|^ = t^- 

As a consequence, for every x G [—s, s], 

/ (^) 


|dz| 


^ 4 
^ 4 


( 3 ) 

rP Afi/^jt— at ^^ .)'^) f4+WA|t|5 

/ e “"JV 4! T M 

/a''-1/4 

rP Ari/4jt—^^^4 


df 


JAr-1/4 

4 ^ 

e dt ^ 


C 


Ari/4 ^ jVi/4 

where the last inequality follows from the fact that liminfArUAr > 0. Similarly, 


(4.27) 




/t+ut- 


(3). S 

^ AT^T J.5 —Af ’^^ ., ^’ F+AfA|t|3 

^ 4iV / Pe ‘^N 4! ^ ' dt 

JAf-1/4 

^ AN Pe ‘^N ^ dt 

yAf-i/4 

= —= / ^ dt ^ —= . 

VN Ji Vn 


(4.28) 


Finally, we consider the integrals over the contours and . By definition 2 : G T+ U 

if and only if there exists 6 G [—7r/4,7r/4] such that 2 ; = cat =t and hence 

|lKe(2-C7v)| ^ |lHe(2-cjv)^| ^ |iHe(2-Civ)^| ^ N-\ |z-Civ|® ^ 
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As a consequence, for every x G [—s, s], 


/ (^) 

e '^iV 2 4! 


/x+uTy 

< g^+^^Jv('^iv)+9^\c]v)+^^^ 


/x+uTy 


|dz| 




c 


Ari/4 ^1/4 , (4.29) 

where for the last estimate, we used the fact that lim inf tv ctat > 0, the decay assumption 
\/^ ffTv(civ) —)• K G M and the convergence g^^\c]\f) —^ Similarly, 


/ (2) 


/ Af|2;-CTv|^e 

'x+uxy 

^ g^+^9Whiv)+5(^’(civ) + ^ ^ C* 


Vn Vn 

By gathering (4.24)-(4.30), we thus obtain the estimate 


|dz| 

(4.30) 


A ^ / |dz| / |drc| e 

X * ^ 




Dn{z,w) 




C 


Combined together with (4.22), this finally yields 


sup 

a;,j/g[—s,s] 


f -I f , 1 -TVi/42,<^_Ovi4_Tvi/4yhGXJv) 

/ dz dw -e <^N Dn{z,w) 

/x* 7s* w' - 




C 

AV4 


Hence, (4.21) is proved, which in turn implies (4.20). Proof of Lemma 4.3 is therefore com¬ 
plete. □ 

The next estimate completes the previous lemma by showing one can replace the constant 
Ng'^{c]sf) in (4.20) by —VNr/ajf, and that the resulting kernel is the Pearcey kernel, up to 
a negligible correction term. 

Lemma 4.4. For every s > 0, 


Ari/4 r 

(2i7r)2f7Ar 7 t, 


f 1 _Tvi/4a,(£^+TVi/4yhGl£M 

dz / dw -e '^N 


w — z 


X giV®(^(z-c^)2+Arfe^(2-c^)4-TV^%^(4«-c^)2-JV®i^^(«;-c^P 


43)/ 


converges as A —)• oo towards Kp^{x,y), uniformly in x,y ^ [—s,s]. 

Proof. Let s > 0 be fixed and recall the definition (2.21) of Kp’g^(x,y). We first show that, 
uniformly in x, y G [—s, s], as A —)■ oo. 


(x, y) = 


dz 


' iN^l^p! C 7 jv 


'sn|hKi^} J-im/*p/a^ 


w — z 


+ 0 


iVl/4 


. (4.31) 
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Indeed, keeping in mind that the sequence {crjy) is bounded, it follows from the defini¬ 
tion (2.19) of S together with simple estimates that for every x,y £ [—s,s]. 




r 1 o 4 o 4 

dz / drc_ 


w — z 




r>-|-CXD 




/ -l-oo 

-oo 


e^l‘l-Vdt -"rdn ^ 


C 


iVl/4 


and 


dz 


^ p—iN^^^p/aj^ pioo \ 

J —ioo JiN^/^o/ajv j 


dw _ - _ 


iN^/^p/aj^ J W Z 




V2 


gn 


iVl/4p 


f-l-oo 4 

I '*'-S-l-r-l- 


/•-l-oo ,4 n 

e^|t|-xdt+-e'^+'+t ) 2 / e-^'^-^du ^ 

vr ) .iN^/^p/a^ iVV4 


from which (4.31) follows. 

As a consequence, by performing the changes of variables ^ i—)• N^^^{z — cn)/o'n and 
w I—>• — cn)/o'n in the right hand side of (4.31) and using the definition (2.25) of a^, 

we obtain 


kS (x, y) = 


(2i7r)2cTAr 7x, 


f 1 _jvi/4a.(£^£i^+jvi/4 

dz / dm-e <^N 


w — z 


o (2) f \ o (3) / \ 


- Pf 

X I e "'"Jv 


+ 0 


iVl/4 


. (4.32) 


In order to complete the proof of the proposition, we now prove the estimate (4.32) still holds 
true after replacement of the constant —^/Nt/g]^ by Ng'j^{cN)- This amounts to showing 
that 


Ari/4 


(3) (’3) 

, /■ , 1 -jvl/4,^(iZ£i^+jv giV )4-^jvl/4 qf!^_jv 9iV ^4 

dz / dm-e ^ 4 ! v «; 


w — z 

¥ 

X e -e 




(4.33) 


converges to zero uniformly in x,?/ G [—s, s], which we now establish by using the same type 
of arguments as in the proof of Lemma 4.3. To do so, we set 

and observe that, since we have the convergences g^^\cN) —>■ ^^^^(c) and '/Ng'^{cN) —)• 
Aat = o(l) as N ^ oo. First, since |z — m|“^ ^ \/2N^/‘^ for every z G T* and m G H*, (4.33) 


An = - 


2g‘Ij 


= 1 

2 2 
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is bounded from above by 


(3) ^3^ 

TT * L^* 




- \/]Vt +VNt 

O.n-^ 


(4.34) 


Next, we use inequality (4.23) with 

N 

„ = jv9kM(^_c„)^-jv»kM(„_c„)y 

SO that \u — v\ ^ '/N/S.i^{\z — cat^ + |t(; — CArp), in order to obtain that (4.34) is bounded 
from above by 




L 


f (3) 

^l/4^Mc(m—OyO _Af ^iVpiv) ^AT ( ^He(w—Cj^)‘^ 


rJN 


|dt(;| 


/ (^) 
4t* 

X [ NAn\w - 


By performing essentially the same estimates as in (4.25)-(4.30), we then prove that (4.34) is 
a O (Atv) as a —)• oo uniformly in x,y £ [—s, s]. This completes the proof of Lemma 4.4. □ 

We have now completed the local analysis around cn. More precisely, by gathering Lem¬ 
mas 4.3 and 4.4, we have established the following result. 

Proposition 7. For every s > 0, 

f f da- ^ ^-n^^^^^-^+n^/V-^+nMz)-nm^ 

{2i'KfaN Jr, Je, w-z 

converges towards K|)g^ (x, y) uniformly in x,y £ [—s, s] as N ^ oo. 


4.3 Step 2: Existence of appropriate contours 

In this section, we prove the existence of appropriate completions for the contours T“ U T“, 
Tp U Tx and H* into closed contours, on which the contribution coming from 
in the kernel KAr(x,y) will bring an exponential decay as N increases to infinity. More 
precisely, we establish the following. 
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Proposition 8. For every p > 0 small enough, there exist N-dependent contours T , T"*" 
and H which satisfy for every N large enough the following properties. 


(1) (a) T encircles the \ - ^’s smaller than cn 

(b) T+ encircles all the \~^ ’s larger than cn 


(c) H encircles all the A ■ ^’s smaller than cat and the origin 


( 2 ) 


(a) T- = T-UT-UT-, 


(b) T+ = T+UT+UT+, 

(c) U ^res 

(3) There exists K > 0 independent of N such that 

(a) lHe(/ 7 v(z) - /at(cat)) ^ -K for all z G T“ ^ U T+3 

(b) 9te(/Ar(w) - /Ar(cAf)) ^ K for all w ^ ^res 

(4) There exists d > 0 independent of N such that 

inf {|z — wl : z € U w £ E] ^ d 


inf {1 2 : — 


w\ 


G T U T'*", w G Hres} ^ d 


(5) (a) The contours T , T"'' and H lie in a bounded subset of C independent of N 

(b) The lengths ofT~, T+ and E are uniformly bounded in N. 

In order to provide a proof for Proposition 8, we use the same approach as in [17, Section 
4.4], from which we borrow a few lemmas. Introduce the asymptotic counterpart of IHc/at, 
namely 

iRtf{z) = —a'iRz{z — c) + log \z\ — J log |1 — xz\ i^(dx), z £ C\IC, (4.35) 

where we recall that /C has been introduced in (4.6). We shall use the following property. 
Lemma 4.5. yiefw converges locally uniformly to Die/ on C\IC, and moreover, 

lim lHe/^(c^) =lHe/(c). (4.36) 

N^oo 

Proof. See [17, Lemma 4.7(a)]. □ 

We now turn to a qualitative analysis for the map 94e/. To do so, introduce the sets 

= {z G C : 94e/(z) < 9Ie/(c)}, 11+ = {2 G C : lHe/(z) > lHe/(c)}. (4.37) 

The next lemma encodes the behaviours of IHc/ as \z\ —)• 00 . 

Lemma 4.6. Both n+ and have a unique unbounded connected component. Moreover, 
given any a £ (0,7r/2), there exists R> 0 large enough such that 


■{z£C: 


TT 

2 


< ^ - a} C , 

(4.38) 

37r 1 


:y-«|cl2+. 

(4.39) 


29 


Proof. See [17, Lemma 4.8]. 

Next, we describe the behaviour of 91e/ in a neighbourhood of c. 
Lemma 4.7. There exist r] > 0 and 0 > 0 small enough such that, if we set 


□ 


Ak = \ z G 


0 < jz — c| < ??, 


vr 


arg( 2 ; — c) — k — 


< 9 


}. 


(4.40) 


then 


^±1) -^±3 C fl-, Ao, A±2, A4 C r 2 +. 

The regions are shown on Figure 6. 

Proof Let 77 > 0 be small enough so that B{c, 2r]) C C \ /C. In particular, one can choose a 
determination of the logarithm such that the map 

f{z) = —a(z — c) + log( 2 ;) ~ 1 j log(l — xz) v{dx) 

is well-defined and holomorphic on B{c, 2r]), and its real part is given by (4.35). Moreover, we 
have f'{z) = g{z) — 0 . Since a = g{c) and g'{c) = g"{c) = 0, a Taylor expansion of / around 
c then yields for every z G B{c,r]), 

2 t\cf{z)-yitfic) - ^g^^\c) 2 t\t{z - c)^ ^ f{z)-f{c)-^g^^\c){z-c)‘^ 

^ wU — cP max 

5 ! ' 


Since iHe(z — c)^ = (—l)^r^ when z = c + because g^^\c) > 0, the lemma follows 

by choosing rj and then 9 small enough. □ 

Let Tl 2 (i+i be the connected component of n_ which contains A 2 £+i. Similarly, let TI 211 be 
the connected component of 0+ which contains ^ 21 - We now prove that the following holds 
true. 


Lemma 4.8. 

(1) We have fli = n_i, the interior of Qi is connected, and for every 0 < a < 7r/2 there 
exists i? > 0 such that 

G C : \z\ > R, —^+a< aTg(z) < ^ — aj C fli . (4-41) 

(2) We have ^2 = ^- 2 , the interior of 0,2 is connected, and there exists R> 0 such that 

G C : jzj > i?, ^ + a < arg(z) < ^ — a| C O 2 . (4.42) 

(3) We have Ila = O- 3 , the interior of Ila is connected and there exists 5 > 0 such that 
B{0,S) C O 3 . 

To prove the lemma, we use the following key properties. 
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Lemma 4.9. (1) If is a connected component of 0+, then ri* is open and, if ri* is 

moreover bounded, there exists x G Supp(i^) such that x~^ G fi*. 

(2) Let 12* he a conneeted component of 12_ such that 12* ^ M. 

(a) If 12* is bounded, then 0 G 12*. 

(b) If 12* is bounded, then its interior is connected. 

(c) //O ^ 12*, then the interior o/12* is connected. 

Proof. See [17, Lemma 4.11], □ 

Proof of Lemma 4-8. We first prove (2). Since 122 is by definition a connected subset of 12+, 
Lemma 4.9-(l) yields that its interior is connected (since 122 is open). Next, we show by 
contradiction that 122 is unbounded. If 122 is bounded, then Lemma 4.9-(l) shows there 
exists X G Supp(i^) such that x~^ G 122 . If < c (resp. x~^ > c), then it follows from the 
symmetry lHe/(z) = 91e/(z) that 122 completely surrounds I 23 (resp. 12i). As a consequence, 
123 ^ M (resp. 12 i ^ M) is a bounded connected component of 12 _ which does not contain 
the origin, and Lemma 4.9-(2a) shows this is impossible. The symmetry IPi.zffz) = D\ef{z) 
moreover provides that 12_2 is also unbounded, and ( 2 ) follows from the inclusion (4.39) and 
the fact that 12+ has a unique unbounded connected component, see Lemma 4.6. 

We now turn to (1). Since 122 is unbounded, then 12i does not contain the origin and it 
follows from Lemma 4.9(2a)-(2c) that 12i is unbounded and has a connected interior. Then, 
(1) follows from symmetry 91e/(z) = 91c/(2), the inclusion (4.38) and the fact that 12_ has a 
unique unbounded connected component (cf. Lemma 4.6). 

Finally, since I 23 is bounded as a byproduct of Lemma 4.8-(2), it has a connected interior 
(Lemma 4.9-(2b)) and contains the origin (Lemma 4.9-(2a)). Moreover, since y{ef{z) —)• —00 
as 2 ; —>■ 0, the origin belongs to its interior and, because of the symmetry d\ef{z) = iRzf{z), 
necessarily I 23 = 12 _ 3 . Hence (3) is established. □ 

Proof of Proposition 8. Given any p > 0 small enough, it follows from the convergence of cat 
to c that for every Nq large enough the points ctvq + po™!^ and cnq + pe~'^'^l^ belong to Ai 
and A_i respectively. Thus both points belong to 12i by Lemma 4.8-(l). As a consequence, 
we can complete the path T+ U into a closed contour with a path T+^(Ao) lying in the 
interior of 12i. Since T^^(Ao) lies in the interior of 12i, the convergence cn ^ c moreover 
yields that we can perform the same construction for all N ^ Nq with T^^(A) in a closed 
tubular neighbourhood T C f2i of T+,,(Ao). By Lemma 4.8-(l) again, we can moreover 
choose T^,,(iVo) in a way that it has finite length and only crosses the real axis at a real 
number lying on the right of 1C. By construction, this yields that the set T is compact and 
that the T+^(A')’s can be chosen with a uniformly bounded length as long as N ^ Nq. Since 
Hi C H_ there exists K > 0 such that iRzf^z) ^ 91e/(c) — 3K on T. Since moreover D^efiy 
uniformly converges to Die/ on T and 91e/7v(cAr) —)• fHe/(c) according to Lemma 4.5, we can 
choose Nq large enough such that iRzfN ^ 51e/ + AT on T and lHe/(c) ^ lHe/Ar(cAr) + K. 
This finally yields that 91e(/Ar(2;) — /Ar(cAr)) ^ —K for all 2 ; G T and proves the existence of 
a contour T"*" satisfying the requirements of Proposition 8 , except for the point (4). 

Similarly, the same conclusion for H follows from the same lines but by using H 2 instead 
of Hi and Lemma 4.8-(2). 

We now turn to the contour T“. Given any p > 0 small enough, for every Nq large enough 
the points cnq —and cnq — belong to A 3 and A _3 respectively. Thus both points 
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belong to ^3 by Lemma 4.8-(3), which contains the origin in its interior. As a consequence, 
we can complete the path T~ U T“ into a closed contour with a path T~^g{NQ) lying in the 
interior of fls, which encircles all the A~^’s smaller than ctv for every N large enough, since 
we assumed lim inf^v-^-oo > 0. Then, we can follow the same construction as we did for T+ 
in order to prove the existence of a contour T“ satisfying the requirements of Proposition 8 , 
except for the point (4). 

Finally, the item (4) of Proposition 8 is clearly satisfied by construction since the sets n_ 
and are disjoint, and the proof of the proposition is therefore complete. 

□ 



Figure 6 ; Preparation of the saddle point analysis. The two dotted paths in correspond 
to and while the dotted path in 14+ corresponds to Eres- 


Equipped with Proposition 8 , we are now in position to establish the remaining estimate 
towards the proof of Theorem 5. 

Proposition 9. For every s > 0, the following quantity 


V/x+UT- is Jt. Je., J W - Z 
converges towards zero, uniformly in x,y G [—s,s], as N ^ 00 . 
Proof. Let s > 0 be fixed. This amounts to showing that 


JV>/4 


dz t du,— 


I- res^ res 


w — z 


and 


iVV4 / dz / dw 


Jr 


w — z 


(4.43) 

(4.44) 


32 






























converge to zero uniformly in x, 7/ G [—s, s]. First, Lemma 4.2 yields 

f gA''(5te/iv(u>)-9tc/]v(cAr))|j^^| ^ f gArg^(cAr)?ie(2-cjv)2+Afg^'(cAr)?ie(2-cjv)"‘+A''L-Civ|® 

7 t , JTt 

and hence, by using the same estimates as in (4.27) an (4.29), we obtain 


,Ar(9rc/^(2)-?Rc/jv(civ)) 


/T* 


|dz| ^ 


Similarly, using instead (4.25), we get 

f g-iV(9te/^W-9tc/jv(cjv))|(|^| ^ 


c 

iVV4 

c 


iVl/4 • 


(4.45) 


(4.46) 


Next, Proposition 8 (5-b) yields the existence of L > 0 independent of N such that, for every 
z G T U H, we have \9\e{z — CAr)| ^ L. Since 

^NfN{z)-NfM{w) ^ gAr(/jv(2)-/iv(cjv))g-7V(/^(t«)-/jv(civ))^ 

by using moreover Proposition 8 (4) and then Proposition 8 (3-a), we readily obtain 


Ari/4 


*- res ^ T-es 


i 


dz i du.— 


W — Z 


€ 


__ e '^N I e 


I 


N{d\tfNiz)-9\cfNicN)) 


1 — I jt 

res^ res 


|d-z| 


,-Af(9rc/iv(«))-*Hc/jv(cjv)) 


|du)| 


jVV4 

^ d ^ 


|dz| 


g-Af(«c/jv(«j)-?Rc/jv(cjv))|^y^| _|_ g 


-TVK 


Idtc 


— I JT^ 

res'---' res 

Because of Proposition 8 (5-b) and (4.46), this yields the exponential decay of (4.43) to zero 
uniform in x, y G [—s,s]. The same holds true for (4.44) by using the same arguments, up 
to the replacement of Proposition 8 (3-a) by (3-b), and hence the proof of Proposition 9 is 
complete. □ 

Finally, we can now easily conclude. 


4.4 Conclusion 

Proof of Theorem 5. Recalling the integral representation (4.4) for Kjv(a:,y), we split the 
contour F into two pieces F“ and F"*", where F“ encircles the A~^’s smaller than cat, and F"*" 
the A~^’s larger than cat. Then, we deform 0 so that it encircles F”. This does not modify 
the value of the kernel. Indeed, there is no pole at w = and the residue picked at w = z 


reads 




2 iTTO']\[ 

and thus vanishes by Cauchy’s theorem since the integrand is analytic. 

Now, we can deform by analyticity T"*", F” and 0 into the respective contours T+, T” 
and H provided by Proposition 8. Hence, KAr(x,7/) equals to 


A71/4 
( 2 z 7 r)^(TAr Jt+uT 


1 _Afl/43, (^-'^iv) _|_Afl/4^ (“-Civ) 


dz (t dw^^e ” ' " -N " -N e^/ivC)-Af/ivM^ 

w — z 


and the uniform convergence (4.3) follows for every s > 0 from Proposition 7 and Proposition 
9. The proof of Theorem 5 is therefore complete. □ 
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5 Hard edge expansion: Proof of Theorem 6 

In this section, we use basic properties of Fredholm determinants and trace class operators. 
We refer the reader to [32, 16] for comprehensive introductions, see also [17, Section 4.2] for 
a quick overview. We also use well-known formulas and identities for the Bessel function that 
can be found in [13]. 

5.1 Preparation and proof of Theorem 6 

Recalling the definition (2.33) of and the definition (4.1) of the kernel Kjv, we consider 
here the scaled kernel 

= ( 5 ' 1 ) 

and denote by Kw the associated integral operator acting on L^(0,s). Thus, if Xmin is the 
smallest random eigenvalue of Mjv (see (2.32)), we have 

^ s) = det {l -K^) 

Recalling the definition (2.29) of the Bessel kernel, [17, Proposition 6.1] states that, for every 
s > 0, we have uniformly in x, y G (0, s) as N ^ oo, 

KNix,y) = K[3"j(x,?/)+ 0 (5.3) 

Notice that the limiting kernel appearing in the right hand side is exactly (2.30) without the 
pre-factor (y/x)"/^, and hence is bounded on (0, s) x (0, s). The integral operator associated 
with this kernel is where E acts on L^(0, s) by E/(x) = x"/^/(x). Erom (5.3), it 

is easy to derive the convergence 

det {I - ^ det (/ - = r„(x) . (5.4) 

Indeed, (5.3) yields det(/ — KAr)L 2 (o,<j) —)■ det(/ — EK^^E“^)£, 2 (o,s) and, for any trace class 
operators A,B on T^(0, s) we have 

det(/ - AB)i 2 (o, 5 ) = det(/ - BA)i 2 (o,s) . (5.5) 

Thus, (5.4) follows because both the operators E and are well-defined on L^(0,s) 

and trace class when a ^ 0, and so are the operators Exj^^ and E“^ when a < 0. See e.g. 
[17, Section 6] for a proof. 

To prove Theorem 6, the first step is to improve (5.3) by making explicit the 1/N- 
correction term, and showing the remainder is of order 1/A^^. More precisely, with defined 
as in (2.33), we prove the following kernel expansion. 

Proposition 10. For every s > 0, we have uniformly in x,y € (0, s) as N ^ oo, 

^N{x,y) 

^ (y) - ^^^{(^MVx)JaiVy) + (x -^ ■ 

(5.6) 
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The proof of the proposition is deferred to Section 5.2. If Q is the integral operator on 
L^(0, s) with kernel 

Q{x,y) = ■^{aJa{Vx)JaiVy) + ix- y)K^^]{x,y)y (5.7) 

then Proposition 10 yields the operator expansion (for the trace class norm on L^(0, s)), 

K«=E(K£>-iQ)E-> + o(T), (5,8) 

Now, to prove Theorem 6, one just has to plug (5.8) into the Fredholm determinant (5.2), 
expand it, and identify the 1/N order term. For the last step, we need the following lemma, 
which essentially relies on a formula established by Tracy and Widom. 

Lemma 5.1. We have the identity 

n- ((/ - K^".>)-'Q) = - (^) si logdet (/ - . 

Proof. Introduce the trace class operator M/(x) = xf{x) on T^(0, s), so that one can write 

Q = ^ {aJaU- )®Ja{V~-) + [M, k£)]) , 

where [A, B] = AB — BA. If R(x,y) is the kernel of the resolvent K^^(/ — then we 

have 

^ logdet (/ - K^Be)LP 0 ,s) = ^)' 

Equations (2.5) and (2.21) of [33] then provide the identity 

4sR(s, s) = (J„(^A), (/ - k£V V„(^A)>^2(o,,) • (5.9) 

Since the right hand side of (5.9) equals Tr ((/ — jQ,(-yA) (g) Ja{V ~'))) the lemma would 

follow provided that 

Tr((I-Kg)-i[M,<j])=0, 

but this is obvious since Tr(AB) = Tr(BA) and commutes with (/ — K^^)“^. □ 

It is now easy to prove Theorem 6. 

Proof of Theorem 6. The identity (5.2) and the operator expansion (5.8) yield 

n*(jV"«7A,i„,„ > s) = det {I - E(k£I - iQ)E-')^,,„„ + o . (5.10) 

By plugging the identity xJ'^{x) = aJa{x) — xJa+i{x) into (2.29), we obtain 

^ _ VxJa+i{Vx)J»{^) - ^J»+i{^)Jc.{Vx) 

^BeWiy) — 
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Combined with the asymptotic behaviour (which follows from the definition (2.28)), 


'^a(\/®) — 


(sign(a))" f ^/x 


a ! 


|a| 


— j (l + 0{x)) , X —> 0+ , 


it is then easy to check from (5.7) that both EQ and QE ^ are trace class and thus, by (5.5), 


det {/ - E(kW - = det {/ - k£' + 

Finally, using the following expansion of a Fredholm determinant: 

det(/ -B) = l- Tr(S) + 0{\\Bf) , 


(5.12) 


where B is trace class with trace norm ||i?|| < 1, we obtain 
det (I - k£> + = det (I - det (I + - Kg)-'Q) 

= det (I - (l + i T,- ((/ - k£>)-1Q) + O (JV-2) ) 

(5.13) 

and Theorem 6 follows by combining (5.10), (5.12), (5.13) together with Lemma 5.1. □ 

We now turn to the proof of the proposition. 


5.2 Proof of the kernel expansion 

Proof of Proposition 10. First, by using the representation (5.11) of the Bessel kernel and 
then the identity 2aJa{x) = xJa+i{x) + xJa-iix), we have 

aJa{Vx)Jaiy/y) + {x - 7/)K^“j (x, ?/) 

= ^(2ajQ,(V^)J«(^) + y/xJa+l{Vx)Ja{\/y) “ ^/V Ja{^/x ) Ja+l{^/y )) 

= ]j^WxJa+l{^/x)Ja{^/y) + ^/yJa{^/x)Ja-l{^/y))■ 

Thus, to prove Proposition 10 is equivalent to show 

KAr(x,y) = K|3“j(x,y) 

- (V^Ja+i(v/^)Ja(Vy) + +0 (^) (5-14) 

uniformly in x, y G (0, s) as N ^ oo. To do so, let 0 < r < i? < liminfjv Ai/2 and introduce 
the map 

Recall that Assumption 2 yields lim inf tv Ai > 0 and thus one can choose a determination of 
the logarithm so that Gn is well-defined and holomorphic on {z G C, |z| < i? -|- 1} for every 
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N large enough. In the proof of [17, Proposition 6.1], the following representation has been 
obtained 


KNix,y) = 


1 


dz 




{2iTT)'^ J\z\=r Z J\^\=R w Z-w\w^ 

Recalling the definitions (2.33) of aM and Cat, straightforward computations yield 

-1 


(5.15) 


G'j^(z) — 


— Y 


N‘^(T]\f KNcTj^f ^ 

J=1 


1 " 


-2 


r) 




Na^ 


-3 


G'NiO) - 


G%iO) = - 


Cn 


8iW2 


N 


and hence we have the Taylor expansion 


Gn{z) — G7v(0) — —^ — 


+ 0 \ — 

AN 16iW^ V 


(5.16) 


uniformly valid for \z\ ^ i? + 1. Plugging this Taylor expansion into (5.15) and recalling the 
contour integral representation (2.30) of the Bessel kernel, we readily obtain as iV —)• oo. 


^N{x,y) = 


1 


dz 


drc 1 


{‘2iTr)^ J\z\=r Z J\^,\=R W Z-W\W 


— ) 6 Z~^ w' A 4 


dz 


drc 1 


(2i7ry Jlzl=r Z w z-w\w 


Z _^\y_\z_w 
— I 6 Z~^ w~^ 4 4 




+ 




Cn \ 1 / dz/ dw 

-^ 7 -777 0 — 0 —(z + rc) — e 4 

161 Vfj^y ( 2 ^ 7 r )2 z J\w\=R w ^ \wJ 




(5.17) 


uniformly for x,y € (0, s). In the light of (5.14), we are left to show that 


( 2 z 7 r )2 J^z:\=r Z J\^\=R W \w/ 


Z\^ _£_l__y _l_2_t£ 

g z~^ w~^ 4 4 


= -2l" 


a/2 


(^^/xJa+li^/x)Jai^/y) + y/yJaiVx)Ja-liVy)^ (5.18) 
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in order to complete the proof of the proposition. But this easily follows from the contour 
integral representations of the Bessel function, 


'^ck(v^ ) - ^ . 


JL — B. , 
e-w 4 cl re 


2-72^/ + = 

2^7^(2^)" z 2z7r ru" w 

(5.19) 

see for instance [17, Eq. (164) and (165)] and perform the respective changes of variables 
z —z~^ and w i—)> —w~^. The proof of the proposition is therefore complete. 

□ 

Remark 5.2. By extending the Taylor expansion (5.16) to higher order terms, the compu¬ 
tation (5.17) easily yields the kernel expansion as N ^ oo, uniform in x,y € (0,s), 


KNix,y) = - 

\y 


a/2 


(x, y) + Y^ -^e^N y )} 


+ 0 


1=1 


NL+i 


for every L ^ 1. The kernels Q^^{x,y) can be expressed in terms of sums of Bessel func¬ 
tions (by using (5.19)). By plugging this formula into the Fredholm determinant (5.2), and 
expanding it, this yields an asymptotic expansion of the form (2.35), although we are not 
able to provide a simple representation for the coefficients cjyl(s) when £ ^ 2. It would be 
interesting to identify these coefficients in terms of Fa(s) if it is possible. 
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